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ABSTRACT 

Lower limits on the power emitted in ultrahigh-energy cosmic rays (UHECRs), which are assumed 
to be protons with energy > 10^^-10^° eV, are derived for TeV blazars with the assumption that 
the observed TeV 7 rays are generated due to interactions of these protons with cosmic microwave 
photons. The hmits depend on the spectrum of the injected UHECR protons. While for a —2.2 
injection spectrum, the lower limits on the powers emitted in UHECRs by lES 0229+200, lES 1101- 
232 and lES 0347-121 are lower than their respective synchrotron luminosities (~ 10'*^ erg s~^); in 
the case of lES 1426-1-428 it exceeds the corresponding synchrotron luminosity by up to an order of 
magnitude. The proposed Auger North Observatory should be able to detect 4 x 10^^ eV cosmic 
ray protons from lES 1426-1-428 within a few years of operation and test the TeV 7-ray production 
model by UHECR energy losses while propagating along the line-of-sight, or constrain the intergalactic 
magnetic field to be larger than ~ 10~^^ G in case of no detection. The lower limits on the apparent- 
isotropic jet power from accelerated lO^^-lO^*^ eV proton spectra in the blazar jet is of the order of 
the Eddington luminosity of a 10^ Mq black hole for a cosmic-ray injection spectrum —2.2 or harder 
for all blazars considered except for lES 1426-1-428. In the case of the latter the apparent-isotropic jet 
power exceeds the Eddington luminosity by an order of magnitude. For an injection spectrum softer 
than —2.2, as is required to fit the observed cosmic-ray data above ~ 10^^-10^^ eV, the Eddington 
luminosity is exceeded by the lower limits on the jet power for all blazars considered. 
Subject headings: cosmic rays - galaxies: active - gamma rays: galaxies - gamma rays: general - 
gamma rays: ISM 



1. INTRODUCTION 

The origin and composition of the highest-ene rgy cos- 
mic rays (CRs) are still unknown (see, e.g., IGaisseri 
|2010( ). They are generally thought to be produced 
in extragalactic sources, due to a lack of power- 
ful Galactic sources capable of producing them and 
due to their near-isotr opic arrival dir e ction distribu- 
tion on large scales (jAbraham et al.l 120091 ). The 
Galactic supernova remnants are leading candidates 
for producing CR s below ~ 10^'' eV (see, e.g., 
iPtuskin et al.l [2"010[ ). Above this energy, denoted as 
ultrahigh energy (UHE), observed CRs may origi- 
nate in acti ve galac tic nuclei (AGNs) with blazars be- 
ing fa v ored (IRachen fc Bicrniann 1993; Bcrcz mskv et al.l 
I2006t IDermer et a" .2009^ .Pe'er et al.. .20091) andTor 
m gamma-ray bursts fGRBs: IWaxmanI 120051: IVietril 
[iMa I Wick et al.|[20p : IRazzaaue et a l."2010^ along with 
their low-luminos ity counterparts fMurasc ct al. 200^ 
IWang et al.l[2007h . 

In the context of the AGN/blazar hypothesis, lack 
of a significant correlation between the known position 
of th e sources and the ar rival directions of the UHE- 
CRs (jAbraham et al.|[2009f) is generally explained by de- 
flections of the charged primaries in the Galactic and 
intergalactic magnetic field (IGMF). The value of the 
IGMF is not known but limits exist, < 10^^ G from 
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cosmological observation s (see, e.g.. IBarrow et all 119971 : 
IKahniashvih et al.l [20091 ) . and > lO'^ MO"^^ G assum- 
ing TeV blazars are persistent sources ([Neronov fc Vovkl 
l2010t iTavecchio et"all I2010b() while > lO-^^-lO"^^ G 
assuming TeV blazars are variable ([Dolag et al.l 120111 : 
IDermer et al.ll201H iTavior et al.llMll) . UHECRs prop- 
agating in voids, with very small value of the IGMF, 
may avoid significant deflection and reach us directly 
from their sources, if nearby. For sources beyond ^ 
100 Mpc, the highest-energy primary CRs cannot reach 
us since they lose energy by interacting with the low- 
energy photons of the Cosmic Microwave Background 
(CMB) and the Extragalactic Backgroun d Light (EBL) 
a phenomenon known as t he GZK effect (lGreisen!'1966t 
iZatsepin fc Kuz'minI Il966t ). The secondary particles 
however, can form an electromagnetic cascade and cas- 
cade radiation, mostly Co mpton-scattercd CMB pho- 
tons (see, e.g., IDai et al.l [20 02: Razzaquc ct al. 200^ 
'Ichik i et al.ll2008t ITavecchio et al.ll2010bl : IDermer et ail 
2011), can reach us as if emitted from the source, if the 
IGMF value is sufficiently small. 

Indeed it has been proposed recently that TeV 7-ray 
emission detetcted from distant blazars and showing no 
significant fiux variatio n is secon dary cascad e emission 
by UHECR protons (Es sev fc Kus cnko 2010; E ssev et al.] 
I20iai2011l) . This proposition opens a possibility to esti- 
mate the UHECR power and subsequently the jet power 
of the TeV blazars, since protons are expected to domi- 
nate the energy budget in such a scenario. In an alter- 
nate scenario, non- variation of TeV flux may arise from 
cascade emission in EBL and CMB from > 10 TeV pho- 
tons, which could be distinguished from p roton-induced 
cascade emission at very high energies ([Murase et al.l 
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120111) . These models, however, cannot explain rapid 
variability, often on hours time scale and over broad fre- 
quency range, observed in the bulk population of blazars 
and thought to be related to the size scale of the emis- 
sion region in the jet. The broadband radio to 7 ray 
spectral energy distribution of TeV blazars consists of 
two peaks: A low-energy peak typically in the UV/X- 
ray and a hig h-energ y peak at ~ 100 GeV (see , 
e.g.. iGiommi fc Padovanl 119941: lAharonian et all 120081: 
iTavecchio et al.l i2010al) . While the low-energy peak is 
generally modeled as synchrotron radiation, the high- 
energy peak is often modeled as inverse Compton (IC) 
emission by the syc hrotron-radiating elect r ons (s ee, e.g., 
iCostamante k Ghis ellini 2002; Bot teheii 120071 ) or by 
hadronic emission (jMiicke et al.ii20Q3l) . Apart from ap- 
parent non-variation of the flux at > TeV in a num- 
ber of high-redshift blazars, models of high-energy 7- 
ray emission from the jet often requires very hard in- 
trinsic spectrum in order to avoid rapid flux reduction 
via 77 pair produ ction with EBL p hotons (see, e.g., 
iFranceschini et all [200l iFinke et"all I2010D . Emission 
from UHECR cascade in the EBL-CMB can potentially 
evade these issues (Essey et al. 2010, 2011). 

In this paper we derive lower limits on the UHECR 
and jet powers of the TeV blazars by assuming that the 
total energy loss rate, integrated over energies > 10^^ eV, 
by UHECR protons in the electromagnetic channels 
(P7CMB e+e~; P7cmb 7r° 77, ir^ ^ e^) and 
after propagating to a redshift z^, when the universe be- 
comes transparent to TeV 7 rays, is equal to the bolomet- 
ric TeV 7-ray luminosity calculated from the observed 
spectra and z^. For the EBL models, not co nstrained by 
the Fe rmi Large Area Telescope (LAT) data ()Abdo et al.l 
[20103), the 77 opacity t-,^ k, 1 at « 0.1. For 
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< z-y cascade emission is not reprocessed efSciently 
and is mostly emitted at very-high energies, far above 
the ~ 0.1 TeV threshold of the air Cherenkov telescopes 
but could potentially b e detected by the p roposed HiS- 
CORE 7-ray detector (jHampf et alj|2011l) which would 
be sensitive in the PeV rage. 

Our limits are based on the assumption that UHECR 
protons do not suffer significant energy losses in the 
blazar jet and escape along the direction within the jet 
opening angle. As mentioned earlier, TeV 7 rays that 
are produced in the jet from UHECR interactions, suf- 
fer attenuation in the EBL. Because of this reduction of 
the 7-ray source flux, the total required UHECR power 
is greater if the observed flux is a combination of the 
source and cascade fluxes than if solely coming from the 
cascade. We do not include deflections of the primary 
proton and cascade e+e" in the IGMF, which reduce 
the flux of the secondary radiation and increase the re- 
quired UHECR and jet po wers (see, e.g.. lRazzaque et aD 
120041: iDermer et all 1201 1[ ). Moreover, large > 10"^ G 
magneti c fleld is expecte d in the filaments in the galaxy 
clusters (|Rvu et al.ll2008t ) which can also reduce the cas- 
cade radiation flux when encountered by the primary 
proton and cascade e"'"e~, thus further increasing the 
required UHECR and jet powers. Our limits, though 
model-dependent, are conservative in this sense. They 
will be weaker if other processes, e.g. leptonic emission 
from the jet and/or from the cascade in EBL-CMB, con- 
tribute significantly to the observed TeV data. Based on 



our estimated lower limits on the UHECR fiuxes from 
the TeV 7-ray blazars, we also calculate the detection 
rate of the primary UHECR protons, which travel along 
the line-of-sight from their sources, by the Pierre Auger 
ObservatoryO 

We calculate the energy losses by UHECR protons in 
the CMB in Sec. 2, and estimate the apparent-isotropic 
luminosity of TeV blazars from spectra measured by air 
Cherenkov telescopes in Sec. 3. We derive our limits 
on the UHECR and jet powers in Sec. 4, and calculate 
UHECR event rate based on the limits in Sec. 5. We 
discuss results and conclude our investigation in Sec. 6. 

2. COSMIC RAY SPECTRA AND ENERGY LOSSES 

We calculate the power-law source spectrum (total 
number of particles generated isotropically from the 
source per unit energy and time intervals) with an index 
—K and between the generation energies Eig-E2g, given 
an apparent-isotropic bolometric luminosity (power) of 
cosmic rays Lqr, as 

Lcr(k - 2)£;,v^ 



dN 

^' dEgdt 



1 - [E^g/E^gY 



The observed energy of a CR on the Earth, however, is 
degraded due to adiabatic or redshift (due to expansion 
of the universe) loss and due to collisional losses with the 
CMB and/or EBL photons. The relationship between 
the generated CR en ergy Epiz) at redshift z and the 
observed energy E is (|Berezinskv et al.ll2006D 



dz' , , 
H{z') 



-b„[(l + z')Eg{z% 



(2) 



where Eg{z') = Eg{z){l + z')/{l + z), H{z') is the co- 
moving Hubble parameter and h^iE) = —dE/dt is the 
energy loss rate of an UHECR at the present epoch 
{z = Q). UHECR protons lose their energy by interact- 
ing dominantly with the CMB photons, creating e+e" 
pairs and tt's. These energy-loss rates and their evolu- 
tion with redshift have been calculated by many authors, 
and they are well understood in the 10^^-10^° eV en- 
ergy r ange (see, e.g..lStanev et al.ll2000l : lBerezinskv et al.l 
[2001 IDermer fc Menonll2'009l) ! 

The flux of cosmic rays on the Earth from a source at 
redshift z is then 



J{E)-- 



q{Eg) dEg 

47rd2(l + z) dE 

LcR(At-2)gr' 
47rd2(l + z)3-'^ 



E- 



where the second line follows with an approximation 
Eig <C E2g in Equation ([T]) and when the lowest-energy 
CRs lose energy due to redshift only, i.e. Eig = Ei{l + z). 
The Jacobian f actor dEg/dE follows fr om Equation ([3]) 
and is given bv iBerezinskv et al.l (|2006f) . The comovmg 
distance dc, which is related to the proper distance dp and 
the luminosity distance as d^ = {l + z)dc = {l + z)'^dp, 

^ |http : //www . auger ■ org/] 



UHECR and jet powers of TeV blazars 
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Fig. 1. — Expected UHECR proton spectra on the Earth from a 
source with increasing redshift from the top to bottom according to 
Equation We assumed an apparent-isotropic UHECR source 
luminosity to be 10^^ erg with a power-law spectrum of index 
K = 2.2 in the generated energy range lO^'^-lO eV. Adiabatic 
(redshift), e~^e~ pair and tt creation losses have been taken into 
account in the numerical calculation. 



is defined as (see, e.g.. lHog3ll999l ) 
c f dz' 



Jo 



(4) 



We assume the standard cosmological parameters: flm = 
0.27, JIa = 0.73 and a Hubble constant Hq — 
71 km s"i Mpc"^ 

Figure 1 shows the UHECR proton spectra as would be 
observed from a source at various redshifts but with the 
same apparent-isotropic Iruninosity Lp = lO''^ erg s~^ 
for 10^'' eV < Eg < 10^° eV and for the spectral index 
of the generated protons to he k = 2.2. Note that the 
sharp cutoff at the highest energies > 4 x 10^^ eV is 
due to TT losses (the GZK effect) and the dip below ^ 
4 X 10 ^^ eV but above ^ 10^ ^ eV is due to e+e^ pair 
losses (jBerezinskv et all 120051 ). These spectral features 
shift to lower energy for high-redshift sources. Adiabatic 
losses dominate at the lowest energy in this plot. 

The total bolometric power lost, due to creating e^e~ 
pairs and tt's, by the UHECR protons, generated at a 
redshift z and measured at another redshift zq, can be 
calculated from Equation ([T]) and Equation 1^ consid- 
ering only coUisional losses and by integrating over the 
generation energy in the range Eig~E2g. The total power 
lost in electromagnetic (EM) channels is smaller when 
considering energy losses to secondary charged particles 
and photons only, which can initiate electromagnetic cas- 
cades by interacting with the CMB and EBL photons. 
Below Eg ^ 4 X 10^^ eV, energy losses by UHECRs are 
mostly in the EM channels. Above Eg ^ 4 x 10^^ eV, 
roughly 1/3 of the bolometric power lost, with ~ 50% in- 
elasticity, goes into producing neutrinos. Subtracting the 
neutrino losses we calculate the ratio of the power lost 
in EM channels, due to propagation, to the generated 
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{k-2)E- 
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1 - iE^g/E2gr-^ Ho (1 + zy 



dEgEg'^ 



^_^, {l + z')bo.^M[{l + z')Eg{z')] 

v/r!„,(i + z')3 + r!A 



Here the lower limit of the redshift integration corre- 
sponds to z^ at which sa 1, since efficient reprocessing 
of the EM energy losses by UHECRs to observed TeV 7 
rays can take place only above this redshift. 

Figure 2 shows the fraction /cr as a function of the 
redshift and for different values of the spectral index 
K. The top and bottom panels correspond to the 10^''- 
4 X 10^^ eV and 10^^-10^° eV energy ranges of the gener- 
ated spectra, respectively. The dotteded curves in each 
panel correspond to the present redshift (zq = 0) in 
Equation ([5]). The solid curves correspond to the red- 
shift zo = z-y = 0.1 at which 77 — > e^e^ pair production 



op acity 



„ 1 at T eV, acc ording to the EBL models 

bv lFinke et all (|2010t) and bv iFranceschini et al.l ([2008), 
and those photons can propagate to us without much at- 
tenuation (Fazio fc S teckcr .1970). Figure 3 shows /cr as 
in Fig. 2 but above 10^* eV for comparison. Note that 
the loss fraction increases in this case because a smaller 
part of the UHECR spectrum is considered here, where 
coUisional losses are the most important. As can be seen, 
the loss fraction crucially depends on the index k of the 
generated CR spectrum and is smaller for softer spectra. 

The Lorentz factor of the e^e~ pairs created from 
a 10^^ eV proton is roughly 7e ~ 10^. The 
Compton-scattered CMB photon energy by the pairs is 
(4/3)72£cMB ~ 900(1 z)(7e/109)2 TeV. Subsequent 
absorption of these photons in the EBL, creating e^e~ 
pairs, and Compton scattering of CMB photons give rise 
to an EM cascade. A similar scenario applies for tt'' — 77 
and TT^ — >■ channels. Compton photons of a given en- 
ergy from the cascade can reach the Earth once their 
opacity falls below unity. Thus substantial reprocessing 
of UHECR power lost in EM channels to observed TeV 
7 rays can only take place at z^ > 0.1. We estimate the 
total EM power loss by UHECRs from TeV 7-ray data 
in the next section. 



3. TEV BLAZARS AND GAMMA RAY LUMINOSITY 

We assume that the TeV 7 rays detected from blazars 
at z > Zj 0.1 are reprocessed cascade emission 
from UHE protons with energy 10^'' < Eg < 10^" eV, 
that are generated at the blazars. TeV emission should 
not be rapidly variable in this case since cascade emis- 
sion takes place over a large size scale. There can 
be some variability at an energy much below 1 TeV, 
due to contamination by the source photons. Cur- 
rently there are four such Te V blazars with adequat e 
data, namely l ES 1426-^428 lA haronian et al.' '20031^ 
lES 0229-1-200 (lAharonian et al. 2007c: Perkins et al] 
l2OT0t) . lES 1 101-232 (lAh aronian et^ [20071 and 
lES 0347-121 (jAharonianet al.M2007bf) . These are all 
BL Lac objects. 

Figure 4 shows the very high-energy 7-ray data col- 
lected from these blazars. None of these sources except 
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Fig. 2. — Total bolometric power loss fraction by UHECR pro- 
tons, generated at various redshifts and with different spectral 
index k, propagting to z = (dotted curves) and to 2 = 0.1 
(solid curves). The universe becomes optically thin (r-,,-, < 1) 
to TeV photons due to e+e" pair creation w ith the EBL pho- 
tons at 2 Si 0.1 (see, e.g., the EBL models bv lFinke et al.l 120101 : 
IFranceschini et al.ll2003) Top panel — Energy losses by e+e" pair 
creation, which dominates below 4x eV, with the CMB pho- 
tons. Bottom panel — Energy losses by e+e" pair and vr creation, 
which dominates above 4 X 10^^ eV, with the CMB photons. 
The lowest UHECR energy is assumed to be 10^'' eV. 

lES 1101-232 has been convincingly detectec0, which is 
sensitive in the 0.1-300 GeV range. However data at 
much lower frequencies, Optical-UV and X-ray, are avail- 
able. Emission at these frequencies is widely believed 
to be synchrotron radiation from non-thermal electrons 
in the blazar jet. The bolometric synchrotron luminos- 
ity igyn, calculated as ten times the specific luminosity 
v Lu at the peak of the s ynchrotron emission as given 
in iTavecchio et al.l (|2010aD , is also listed in in Table 1 for 
each of the four blazars. 

Figure 4 also shows t he de absorbed 7-ray data points 
using the iFinke et al.) (|2010D EBL model at = 0.1. 
We fit the deabsorbed 7-ray spectra (black dotted lines) 
by two smoothly-joined power laws, similar to a Band 

TABLE 1 

Blazar and TeV 7-ray model parameters 
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(erg 


h) 




lES 1426 -f 


428 


0.129 


6 X 10''^ 


2.1 
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10^5 


-2.3 


lES 0229 -1- 


200 


0.139 


2 X lO''® 


4.6 


X 


10** 


-1.7 


lES 1101 - 


232 


0.186 


3 X lO''® 


2.1 
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10** 


-2.0 


lES 0347 - 


121 


0.188 


2 X 10*6 


2.0 
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10** 


-2.2 



^ A tentative d etection of lES 022 9 -1-200 has been reported 
bv lUrr et aEI ll20Tll'l and bv lAbdo et"all Il2010bl') 
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Fig. 3. — The same as in Fig. 2 but the lowest UHECR energy 
is assumed to be 10^* eV. 



spectrum (jBand et al.l Il993( ). with a super-exponential 
cutofi" at the high energies as below 



dN 
dedAdt 



= C 




(2+a)£ i . (a-ff)£pk.GcV 

Spk J ' ^ 2+a 

(a-)3)epk,GoV ia~3 (R\ 

2+a J ^ 

X f./3-"p-(e/ec)'' . p > (a-ff)£pk.GcV 
AC e , t ^ 2-i-a 

We keep the index a of the low-energy power law fixed 
such that vFi, o^v^ (a = 0), the hardest Compton spec- 
trum expected, to be conservative. The break energy Cpk, 
at which the vF^, spectrum turns over, the high-energy 
power-law index /3 and the cutoff energy were allowed 
to vary. 

Finally we calculate a conservative lower limit on the 
apparent-isotropic TeV 7-ray luminosity from our fitted 
spectra (see Fig. 4, dashed lines) to the deabsorbed data 
points at = 0.1 as 



iz^=o.i =47rd|(z^) / dv Fj 



(7) 



in the 10 GeV-20 TeV range. This luminosity and the 
high-energy photon index /3 in Equation ([6|) along with 
the synchrotron luminosity for individual blazars are 
listed in Table 1. 

4. LIMITS ON UHECR AND JET POWERS 

We calculate the lower limits on the apparent-isotropic 
UHECR power of the TeV blazars from 7-ray luminosity 
calculated in Equation ([7]) as 



L 



UHECR 



> iz-,=0.l/ /CR 



(8) 



Note that for high EBL models (see, e.g., iStecker et al.l 
|2006() . with z.^ < 0.1, will be lower. A lower Zj 
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Fig. 4. — TeV 7-ray data (filled points) from the HESS and VERITAS observations on selected blazars at z > 0.1. F its (dashed l i nes) t o 
the 7-ray data from these blazars are performed after deabsorption (empty points) at z = 0.1 with the EBL model bv lFinke et al.l II2010I ). 
The fit function consists of two smoothly- joined power laws with a super-exponential, exp(— v^), cutoff at high energies as in Equation ||6ll. 
The low-energy slope of the power law is kept fixed so that uFi, oc u^, analogous to the hardest IC spectrum. The high-energy index f) is 
listed in Table 1. 



also gives a higher /cr or more efficient conversion of 
UHECR energy to EM energy. Thus our hmits wiU be 
weak er for high EBL m odels, a result consistent with 
what lEssev et al.l (|2011 f) found by modeUng TeV 7-ray 
data. 

The top panel of Fig. 5 shows the lower hmits on the 
apparent-isotropic UHECR power, in the 10^''-10^° eV 
range, of individual blazars that we have derived from 
the TeV 7-ray luminosity Lz^=o.i listed in Table 1 
and Equation ([5]). The limits are in the range 10''^- 
lO^'' erg/s and are consistent with the required UHECR 
sourc e luminosity frorn detailed modeling of TeV 7-ray 
data (lEssev et al.ll2010f) . 

The bottom panel of Fig. 5 shows the lower limits on 
the blazar jet power that we have derived from the TeV 
7-ray data, and are plotted against the synchrotron lu- 
minosities of the corresponding blazars. To calculate 
these lower limits we have extrapolated the proton spec- 
trum down to a minimum energy of 10 GeV. Thus the 
bulk Lorcntz factor of the relativistic TeV blazar jets 
need to be ~ 10, which is consistent with measured val- 
ues. The empty symbols (as in the top panel) corre- 
spond to the limits assuming the same power-law in- 
dex K as above 10^^ eV and as is expected in a Fermi 
shock-acceleration scenario. The filled symbols corre- 
spond to the limits assuming k — 2.2 below 10^'' eV 
and K — 2.5 or 2.7 above 10^^ eV. Such broken power- 
law spectra may arise in more complicated scenarios, e.g. 



from cooling of UHECR protons in the jet. In case the 
minimum proton energy is 100 GeV, our limits will be 
weaker by factors 1.6, 3.2 and 5.0 for k — 2.2, 2.5 and 
2.7 respectively, in case of single power-law spectrum. 
Also shown (dotted line) is the Eddington luminosity 
LEdd = 1.3 X 10^^ (Mbh/10^ Mq) erg/s of a 10^ Mg 
blackhole, which is the typical black-hole mass thought 
to power the blazar jets. 

In case the TeV blazars produce UHECRs above 
10^® eV rather than above 10^^ eV, our limits on the 
UHECR power becomes lower (see Fig. 6 top panel) by 
a factor of ~ 2-5, depending on the injection spectral in- 
dex K. This is an expected result since the blazars need 
to power a smaller part of the UHECR spectrum in this 
case. The lower limits on the jet power (Fig. 6 bottom 
panel) in case of broken power-law spectra (filled sym- 
bols) are also reduce by a factor ~ 2-3 while shifting the 
break energy from 10^^ eV to 10^® eV, for similar reasons 
as for the limits on the UHECR power. Our limits on the 
jet power for single power-law spectra (empty symbols in 
Fig. 6 bottom panel) from 10^" eV to 10^° eV, however, 
remains unchanged. 

The limits that we have derived are based on the as- 
sumption that primary UHECRs and the cascade emis- 
sion are concentrated within the opening angle of the 
blazar jet, 6'jct ~ 0.1. For 10^'' eV CRs and a char- 
acteristic propagation distance of dp « 500 Mpc with 



Ir 



1 Mpc coherence length scale of IGMF, this im- 
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plies that the mean IGMF strength should be 



B 



IGMF. UHECR 



< 10-12 



Oj^\ ( E 
0.1 / V 1017 eV 



500 Mpc 



-1/2 



^coh 

Mpc 



-1/2 



G, (9) 



for the deflection angle to be smaller than ^jot- A stronger 
constraint on the IGMF, however, comes from the re- 
quirement that the cascade electrons of Lorentz factor 
7e « 10^ that upscatter CMB photons to TeV 7 rays do 
not deviate out of The cooling length scale of elec- 
trons in nearby universe is Aic = 0.71 (7e/10^)~i Mpc, 
and to avoid any deflection larger than 0jot one requires 



B 



IGMF.TgV 



< 10-16 ( 11'=* 

0, 



- f^y 

1 / VTeV/ 



G. 



(10) 



To avoid significant deflection in the Mpc scale 
struct ured region of g alaxy clusters with ~ 10~^ G 
field (iRvu et all [20081 iMurase et all 1201 ID . the coher- 
ence length scale should be < 10 kpc from Equation ©. 
For much larger magnetic field and/or coherence length 
the UHECRs become isotropic and their luminosity de- 
creases by a beaming factor /f, = (1 — cos0jot)~i ~ 200 
for 0jet = 0.1. As a result our lower limits on UHECR 
power will become stronger in such a case. 

Propagation of UHECRs from the TeV blazars nearly 
along our line-of-sight, in low IGMF environment, brings 
up naturally the question of their detectability by the 
UHECR detectors on the Earth. While the highest- 
energy CRs are absorbed during propagation in the 
CMB, since the blazars are beyond the GZK radius, suf- 
ficiently energetic CRs can still reach us (see Fig. 1) and 
we discuss the prospect of their detection in the next 
section. 

5. DETECTABILITY OF UHECRS FROM TEV BLAZARS 

Th e location of the Pierre Auger Observa- 
tory () Abraham et al.l 120041 ) at latitude 35.2° S makes it 
insensitive to the nearest and most powerful UHECR 
candidate source lES 1426-1-428 with a Declination 
5 = -1-42.67°. The proposed Auger North detectoiQ 
at latitude 39° N, however, will be sensitive to lES 
1426-F428. For a detector similar to the Auger South, 
the instrument exposure becomes independent of energy 
at > 3 X IQi'^ eV and d epends only on the zenith angle 
H Abraham et all [20081 ) ■ For Auger South 6^^ < 60°. 
The relative exposure at a point source in the sky, such 
as lES 1 426-1-428, comp ared to the largest exposure on 
the sky (jSommerj [200l[ ) is uj{5) = 0.62 for the Auger 
North. The actual exposure on the source can be found 
from multiplying uj{S) by the total integrated exposure 
£ (km^ yr sr) over the detector's field-of-view and 
dividing by the solid angle of the detector O = y/Sn sr, 
for 9^ < 60° (see also lCuoco fc Hannestad[[2008[ ). Thus 
the total number of UHECR events expected from 
lES 1426-1-428, for k — 2.2 source spectrum, in the 
proposed Auger North detector is 



n 



J(E) dE 



Eth 



' |http : //www ■ augernorth . org/ 1 



1046 erg/s J \9, 000 km2 yr sr^ 



(11) 



above a threshold energy E'th = 4x 10 eV and below the 
maximum energy E^^ = 10^° eV. Thus lES 1426-1-428 
should be clearly visible, if an UHECR source as modeled 
here, with the same exposure used in the first data releas e 
by the Auger South Observatory ([Abraham et al.|[2008f) . 

The arriving UHECRs from lES 1426-^428 are ex- 
pect ed to be defi ected by the Galactic magnetic field (see, 
e.g.. iDermer et a l. 2009) by an angle 



1.5° 



E 



4 X 1019 eV 



B, 



Gal 



10-6 G 



\ kpc 



sin b 
sin 64.9° 



(12) 



where /idisc ^ 1 kpc is the height of the Galactic disc 
and b = 64.9° is the galacic latitude of lES 1426-^428. 
The deflection angle while propagating in the IGMF, 
with -BjGMF.TcV < 10-16 G, is negligible. Thus UHE- 
CRs from lES 1426-1-428 should be clustered around 
the source direction for a Northern site detector simi- 
lar to Auge r South with ~ 1° or better resolution above 
- 1019 eV (lAbraham et al.l[200l) . 

It is interesting to note that the position of 
lES 1426-1-428 is within - 5°-10° of three clustered 
(triplet) events detected by the AGASA experiment 
during its operation with energie s 4.97 x lOi^ eV , 
4.98 X 1019 eV and 5.27 x lOi^ eV ([Takeda et al.l[l999l : 
iHavashida eTall [20M 1. The average angular resolu- 
tion of AGASA however is 1.8°, smaller than the an- 
gular separation between the triplet and lES 142 6-1-428. 
Down grading the event energies by 25% (Taked cTet al.l 
[2001 does not increase the deflection angle ^Qai sig- 
nificantly. If, however, i?Gai ~ 4 x 10~6 G in the 
Galactic disc, 6'Gai can be large enough to encompass 
the nearest AGASA event from the triplet. Although 
the exposure is lower compared to its peak exposure in 
other directions in the sky, HiRes experiment has not de- 
tected any even t above 4 x IQi^ eV from the direction of 
lES 1426-f428 ([Abbasi et al.l[20Tol ). 

The detection prospect for UHECRs with k = 2.2 from 
the other 3 blazars within the field-of-view of the Auger 
South is much smaller, iVovt = 0.04, 0.03 and 0.02 re- 
spectively for lES 0229-1-200 {6 = -f20.27°), lES 1101- 
232 {6 = -23.50°) and lES 0347-121 (<5 -11.98°); for 
Lp ^ 10"*^ erg/s and £ = 9, 000 km^ yr sr. 

6. DISCUSSION AND CONCLUSIONS 

The lower limits on the UHECR power -LuhecRj in 
the lOi^-lO^'' eV range, that we have derived are lower 
than the synchrotron luminosity Lgyn of the sources for 
K = 2.2, except for lES 1426-1-428 where LuHECR/^syn ~ 
6, 15, and 30 for k = 2.2, 2.5, and 2.7 respectively. While 
none of these sources hav e been detected by the Fermi 
LAT (lAbdo et al.ll2010aD . Luhecr is a factor ~ 1-60 
times the median 7-ray luminosity of ^ 3 x 10'*^ erg/s, 
in the LAT range, from th e 118 BL Lacs de tected by 
Fermi with known redshift ([Abdo et al.ll2010£j) . The ra- 
tio iuHECR/-^syn gets Smaller in the case TeV blazars 
generate cosmic rays in the lOi^-lO^^ eV range only. 
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The lower limits on the jet power that we have calcu- 
lated, exceed Lgyn and the median BL Lac 7-ray lumi- 
nosity in the Fermi LAT range for all the blazars we con- 
sidered. In the most extreme case of lES 1426-1-428 the 
ratio ijct/isyn ~ 2 X 10^, 5 x 10^, 2 x 10^ for the single 
power-law spectrum with n = 2.2, 2.5, 2.7 respectively. 
Even for k — 2.2, the limiting Eddington luminosity is 
violated in this case by an order of magnitude. For other 
blazars, Ljot ^ -^Edd for k = 2.2 and single power-law 
spectrum. However Lj^t ^ L^dd for k > 2.2, both for the 
single power-law and broken power-law spectra. 

If the observed TeV 7 rays from the blazars that we 
have considered are due to UHECR acceleration in these 
sources and due to their energy losses while propagating 
along the line-of-sight, the source spectrum of UHECRs 
needs to be oc E or harder in order for the model 
to be energetically viable. If the synchrotron luminos- 
ity is a good indicator of the contribution by the rel- 
ativistic electrons to the jet power, then the power in 
baryons with Lorentz factor 7p > 1 in the relativistic 
jet needs to be larger by a factor of 200 or more in the 
case of lES 1426-f428 with k = 2.2. The actual jet 
power will be even higher when including non-relativistic 
baryons in the jet. The proposed Auger North Observa- 
tory should be able to detect CRs above 4 x 10^^ eV 
from lES 1426-f428 within a few years of operation and 
test the mode l of TeV 7 ray origin from UHECRs as pro- 
posed by Esse v et all 0010, 2011) and as discussed here. 
Detection of GZK neutrinos from lES 1426-1-428 by the 
currently operating IceCube neutrino observatory at the 
South Pole and by its proposed extensions could be pos- 
sible. Non-detection of UHECRs from lES 1426-^428 
can imply a large IGMF as generally inferred from a 
lack of significant correlation between the sources and 
UHECR arrival directions, and/or favor leptonic produc- 
tion mechanism of TeV 7 rays. 

To fit the observed UHECR spectrum as measured 
by the Auger and HiRes arrays, with protons, an in- 
jection spectrum cx above 10^^ eV is required to 
explain the dip and cutoff as due to e^e~ pair produc - 
tion and tt production, respectively (jBerezins kvl [20081 ). 
A broken power law cx above and oc E~'^ '^ be- 
low a break-energy ~ 10^^-10^^ eV can help reduce 
the energy crisis in the AGN interpretation of the ob- 
served UHECRs. This scenario also makes the en- 
ergy requirement for the UHECR interpretation of TeV 
7 rays from blazars less severe. To match the ob- 
served spectrum, the total energy injection rate in UHE- 
CRs above lO'^^ eV in the local universe needs to be 



^ 10"^^ erg Mpc~^ y r~^ (see e.g. iVietri et al.l 120031: 
iBerezinskv et al.|[2006t) . The energy injection rates above 
10^^ eV from the 4 TeV blazars that we considered are 
2,> 6, > 11) X lO'^Vj^^ erg Mpc'^ yr"!, respec- 
tively, for K, = (—2.2, —2.5, —2.7) from our lower lim- 
its on UHECR power. The energy injection rates are 
calculated from the UHECR powers divided by the co- 
moving v2luines_corresponding to the redshifts of each 
blazar (jPermer fc R azzaauc 2010). If the structured re- 
gion with high magnetic field near a TeV blazar makes 
> lO^'' eV cosmic rays isotropic, our lower limits on 
UHECR power will be higher by a factor ~ fb and the 
corresponding energy injection rate could be sufficient to 
explain the observed UHECR spectrum above ~ 10^^ eV. 
There will be many TeV 7-ray sources without the 
UV/X-ray synchrotron counterparts from the misaligned 
blazars in this scenario, which could be searched for with 
existing and future telescopes. 

The jet power exceeds the Eddington luminosity in al- 
most all cases we considered and will do so more severely 
if the UHECRs become isotropic in the scenario discussed 
above. It is possible that extragalactic sources with hard 
injection spectrum contribute to the observed UHECRs 
far above ^ 10^^ eV, as explored in the GRB model with 
K ^ 2.2 above the ankle at ~ 4 x 10^^ eV and wit h 
observed star-formation rate (jPermer fc Atovan|[200l . 
However, in such models an additional component of 
UHECRs in the intermediate - IQ^'^-A x 10^* eV en- 
ergy range is required. In principle they can originate 
from low-luminosity GRBs/hypernova e of either extra- 
galactic jWa ng et al. 20 03) or Galactic (I Wick et al.ll2004l: 
IBudnik et al.ll2008t iCal-rez et al.l 120101 ) origin. Whether 
nature conspires to generate such a high-energy break 
in the injection spectrum or the observed UHECRs orig- 
inate in multiple source classes is a matter of ongoing 
investigations. However, the interpretation of TeV 7 
rays from blazars generating UHECRs seems less favor- 
able in case the injection spectrum is similar to the one 
required to explain the observed cosmic-ray spectrum 
above ~ lO^^- lO^* eV. 
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